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Abstract 

tJ" . A variant for tlie Hilbert and Polya spectral interpretation of the Riemann zeta 

I function is proposed. Instead of looking for a self-adjoint linear operator H, whose 

' spectrum coincides with the Riemann zeta zeros, we look for the complex poles of 

Q . the S matrix that are mapped into the critical line in coincidence with the nontrivial 

D I Riemann zeroes. The associated quantum system, an infinity of "virtual resonances" 

Q I described by the corresponding S matrix poles, can be interpreted as the quantum 

OO ' vacuum. The distribution of energy levels differences associated to these resonances 

shows the same characteristic features of random matrix theory. 
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1 Introduction 

There is a conjecture, quoted to be made by Hilbert and Polya, that the zeros of 
the Riemann zeta function (^{z) on the critical Hne are eigenvalues of a "mysterious" 
complex Hermitian operator H. We report in this presentation a variant of the above 
conjecture where we associate the unknown quantum system to the "vacuum". This 
vacuum is interpreted as an infinity of "virtual resonances", described by complex 
poles of the scattering S matrix. 

Actually the Riemann zeta zeros fluctuations are interesting by their universal- 
ity, being observed in quantal spectra of different physical systems (see Refs.P), 
and by the connection they have with chaotic dynamics, (for a review see Bohigas 
12). Ever since Montgomery's [Sj conjectured these zeros behave like the eigen- 
values of a random hermitian matrix, attempts have been made in order to find a 
quantum system with the Hamiltonian represented by such an operator (see Berry 
and Keating 4 , and references therein). Our goal here is to present a new conjec- 
ture in order to relate the origin of this universality to the statistical fluctuation of 
the vacuum that provides the same bakground for the different physics systems in 
interaction with it. 

We will begin by a short outline of the concepts of the Riemann zeta function, 
prime numbers and the Jost function. After, motivated by recent work [S] in which 
was shown numerically a one to one connection between the large zeros of Jost 
function in the complex momenta plane with large prime numbers and large complex 
Riemann zeta zeros, we propose a distribution for these Jost zeros representing the 
quantum vacuum. Finally, rather than look for a potential that put these Jost zeros 
on the desired place, we show that the corresponding resonances energies has the 
same statistical fluctuation given by the random matrix theory (RMT). We conclude 
suggesting an approximate model for the above quantum vacuum. 

2 Prime numbers, Riemann's zeta function and 
Jost function 

A comparison between prime numbers, Riemann's zeta function and Jost function 
has been given in a recent paper 0. In this section we summarize it in order to 
introduce the proposed map © between complex zeros of Jost function with the 
complex zeros of the zeta function. For a recent review concerning the relevance of 
prime numbers to Physics see Rosu 

The approximate number of primes 7r(x) less than a given x, also called the 
prime counting function, is given by the prime number theorem tt{x) ~ log x (see 
Titchmarsh [7], Chapter- III). This relation gives the asymptotic approximation for 
nth prime Pn, 

Pn = nlog n , as n ^ oo . (1) 

The connection between the distribution of prime numbers irlx) and the complex 
zeros of the zeta function started with Riemann's 1859 paper (see Edwards jHI, P- 
299). Riemann's zeta function is defined (Ref.[7], p.l), either by the Dirichlet series 
or by the Euler product 

Ciz) = Y.n-^ = n(l - p-')-' , Rez>l, (2) 

n p 
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where n runs through all integers and p runs over all primes. This function can 
be analytically continued to the whole complex plane, except at z = 1 where 
it has a simple pole with residue 1. It satisfies the functional equation (^(z) = 
sm{TTz/2)T{l — z)({l — z) , called the reflection formula, where T{z) is 
the Euler gamma function. It is known that C,{z) has simple zeros at points 
z = —2n, n = 1,2,..., which are called trivial zeros, and an infinity of complex 
zeros lying in the strip < Re z < 1. Due to the reflection formula they are sym- 
metrically situated with respect to the axis Re z = 1/2, and since C{z*) = C*{z), 
they are also symmetric about the real axis. Consequently it suffices to consider the 
zeros in the upper half of strip 1/2 < Re z < 1. It is possible to enumerate these 
complex zeros as Zn = bn + i tn , with ti < t2 < < . . ., and the following result 
can be proven (see Titchmarsh j3| , p. 214) 

27rn 

Fn ~ in ~ 1 , as n ^ oo. (3) 

log n 

The Riemann Hypothesis (RH) is the conjecture, not yet proven, that all complex 
zeros of ({z) lie on the axis Re z = 1/2, called the "critical line". This is usually 
considered as the most important unsolved problem in Mathematics. Based on this 
conjecture, Riemann derived an exact formula for 7r(x) taking into account local 
prime fluctuations in terms of nontrivial zeta zeros (Ref.|H], p. 299). 

On the other hand, the Jost function has played a central role in the develop- 
ment of the analytic properties of the scattering amplitudes. In order to recall its 
properties, let us consider the scattering of a non-relativistic particle, without spin, 
of mass m by a spherically symmetric local potential, V{r), everywhere finite, be- 
having at infinity as V{r) = 0{r~^~'^), e > 0. The Jost functions f±{k) are defined 
(see Newton 0, p. 341) as the Wronskian W, f±{k) = W[f±{k,r), (^(fc,r)], where 
Lp{k, r) is the regular solution of the radial Schrodinger equation 

[^ + k'-V{r)]^{k,r)=0, (4) 

(in units for which h = 2m = 1) k being the wave number and the Jost solutions, 
f±{k, r), are two linearly independent solutions of eq.Q). They satisfy the boundary 
conditions, lim,.^oo[e^*'^^/±(^i ''')] = li corresponding to incoming and outgoing 
waves of unit amplitude. 

The properties of the solution of the differential equation @ define the domain 
of analyticity of the Jost functions f±{k) on the complex fc-plane as well as its 
symmetry properties, such as f^{k*) = f-{k). The phase of the Jost function is 
just minus the scattering phase shift 6{k), that is f±{k) = |/±(A:)| e^*'^^'^) , so that 
the usual S matrix is given by 

Sik) ^ e'^'C'^ = 7=7^- (5) 

The complex poles {Re k ^ 0) of S{k), or zeros of f+{k), correspond to the solutions 
of the Schrodinger equation with purely outgoing, or incoming, wave boundary 
conditions. Resonances show up as complex poles with negative imaginary parts, 
their complex energies being 

kn = — i (6) 
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where H„ and represents the energy and the width, respectively, associated with 
nth resonance state in order of distance from origin. For smah F„, resonances 
appear as long-hved quasistationary states populated in the scattering process. If 
the width is sufficiently broad no resonance effect will be observed and we will call 
this kind of S matrix pole as "virtual resonances" throughout. Resonances, in fact, 
are represented by pairs of symmetrical S matrix poles in the complex k-plane, 
a capture state pole in the third quadrant and the decaying state in the fourth 
quadrant, since they give to the asymptotic solution an incoming growing wave and 
an outgoing decaying wave, respectively, exponential in time |1U| . This could be 
described by Gamow vectors treated by Bohm and Gadella as pairs of S matrix 
poles corresponding to decay and growth states, where the traditional Hilbert space 
description is replaced by the generalized Rigged Hilbert Space in order to account 
for time asymmetry of a resonant process. The transient states which last only 
a very short time corresponding to the above "virtual resonances", to which our 
instruments are insensitive, will be called here as "vacuum". Examples of broad 
resonances are the well known large poles, related basically to the cutoff in the 
potential without any physical interpretation (see Nussensveig p. 178). 

In truncating the potential, i.e., the potential is set equal to zero for r > R > 0, 
which is the cutoff of the potential at arbitrarily large distances R, it is possible 
to obtain explicit asymptotic formulas for the Jost zeros. With this restriction 
it can be shown that equation f-\-{k) = is entire of order 1/2 and according to 
Piccard's theorem has infinitely many roots for arbitrary values of the potential. The 
asymptotic expansion of fc^, for large n, after introducing dimensionless parameter 
13 = kR, is given by (Ref.^j, p.362) 

/3„ = nvr-.i^±4^+0(l) (7) 

where n = ±1, ±2, ±3 • • • and a is defined by the first term of the potential 
asymptotic expansion, near r = R, through V{r) = C{R — r)" + • • •, cr > and 
r <R. 

The connection between the complex zeros of the Jost function and those one of 
the Riemann zeta function is provided by the transformation, 

by which the lower half of complex /3-plane {Re /3 7^ 0) is mapped onto the critical 
axis, Re z = 1/2, of complex z-plane. 

Now we show that for cutoff potentials the transformation (jH)) gives rise to com- 
plex Jost zeros with the same asymptotic behavior as the complex Riemann zeta ze- 
ros, being all in the critical line. After introducing dimensionless quantities, energy 
En = R^'^n and widths G„ = R^T^ , the equation @ is written as = En — iGn/2, 
then by (jEJ we get 

1 En I % 

from ((T)), with a = 0, we see that {Im Zn} has the same asymptotic expansion as 
given by ®, i.e.. 
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which means that for each resonance, the ratio between the energy and the width 
is given by the height of the zeta zero on the critical hne. On the other hand, 
dimentionless widths defined in © as 

G„ = 4 i?e I3n Im , (11) 

after taking into account O, when £7 = 0, shows the same asymptotic expansion 
for large primes given by the prime number theorem, 

Gn = 47rp„ as n ^ oo. (12) 

Then nth large complex Jost zeros are also related to nth large primes, showing 
an asymptotic connection between primes and complex Riemann zeta zeros, in a 
one-to-one correspondence. 

3 Hilbert-Polya conjecture 

The Hilbert-Polya conjecture is the spectral interpretation of the complex zeros 
of the Riemann zeta function as eigenvalues of a self-adjoint linear operator H in 
some Hilbert space. Such H could prove the RH. We suggest a variant for the above 
conjecture: Instead of looking for H, whose spectrum have to coincide with the 
Riemann zeta zeros, we are looking for a potential that gives a Jost function with 
all zeros on the lower half of complex /3-plane {Re /? 7^ 0) that coincide with the 
complex zeros of the Riemannn zeta function after the transformation Q. In this 
way the Riemann hypothesis follows. For real potentials, these complex /3 zeros are 
located symmetrically about the imaginary axis, then by (jSJ they will be mapped 
symmetrically about the real axis into the critical line. The associated quantum 
system could be identified with the quantum vacuum, interpreted as an infinity of 
"virtual resonances", described by corresponding S matrix poles. 

The first to associate the Riemann Hypothesis with transient states were Pavlov 
and Faddeev [T3|, by relating the nontrivial zeros of the zeta function to the com- 
plex poles of the scattering matrix of a particle on a surface of negative curvature. 
Khuri ini has recently proposed a modification to the inverse scattering problem in 
order to obtain the potential whose coupling constant spectrum coincides with the 
Riemann zeta zeros. 

Assuming that the RH is true, we conjecture a vacuum spectrum with widths 
given by G„ = Pn , according to l\V2\i . and the corresponding energies given by 
En = Gntn = Pntm according to (|1(H) . Before trying to find a potential that will 
give these S matrix poles we will study numerically the statistics fluctuation of these 
"virtual levels" . 

4 RMT and energy distribution of virtual res- 
onances 

The random matrix theory (RMT) describes fluctuation properties of quantal 
spectra. In lack of analytical or numerical methods to obtain the spectra of com- 
plicated Hamiltonians, Wigner and Dyson (see Refs. pQ) analyzed ensembles of 



The Riemann Zeta Function and Vacuum Spectrum 



5 



n = 1 ... 72000 



n = 2000 ... 4000 



0.9 - \ 
0.8 - 



n = 10000 ... 12000 




- Vacuum 
GOE 

- GUE 

- - Poisson 




Vacuum 
GOE 
GUE 
Poisson 



n = 70000 ... 72000 




- Vacuum 
GOE 

- GUE 

- - Poisson 




- Vacuum 
GOE 

- GUE 

- - Poisson 



Figure 1: The statistical of level-spacings in vacuum spectra {pnin}- The histogram at 
the top left gives a sequence comprising the first 72000 vacuum levels. The histogram 
at the top right comprises 2000 consecutive levels in the analysis corresponding from the 
2000th to the 4000th level. The histograms at the bottom give 2000 consecutive levels 
centered at different E values. 
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random matrices, in which the matrix elements are considered to be independent 
random variables, that have in common only symmetry properties like hermiticity 
and time-reversal invariance. The RMT predicts, for the nearest-neighbour spacings 
of the energy levels distribution P{s), the form which is described by the Wigner 
surmise: P{s) = A exp{—Bs'^), where Sn = {En+i — En) are spaces between 
adjacent levels of the unfolded spectrum with mean distance < s„ >= 1. It is 
obtained by accumulating the number of spacings that lying between (s, s + As) 
and then normalizing P{s) to unity. A and B are normalizing constants and D 
is a parameter which depends on the symmetry of the system which characterizes 
the repulsion between neighbor levels. For hermitian time-reversal Hamiltonians 
(Gaussian Orthogonal Ensemble) D = I and Pgoe{s) = (vr/2) s exp(— 7rs^/4). For 
complex hermitian matrices, when time-reversal invariance does not hold (Gaus- 
sian Unitary Ensemble) D = 2 and Pgue{s) = (32/7r^) exp(— 4s^/7r). If the 
eigenvalues of a system are completely uncorrelated (Poisson Ensemble ) we have 
Ppe{s) = exp{-s). 

The nearest-neighbor spacing distribution P{s) corresponding to the vacuum 
spectrum energy, i.e. resonances energies given by {pntn}-, where prime sequence 
{pn\ are taken from the table ^Hl and {tn\ computed by Odlyzko J7j are shown 
in Figure 1. The curves corresponding to Poisson distribution (dashed line), GOE 
distribution (dotted line) and GUE distribution (dotted-dashed line) are drawn for 
comparison. For a sequence comprising the first 72 000 vacuum levels we obtain, 
see the top left of Figure 1, a mixing between Poisson and GOE distribution. In 
order to compare with spectra of atomic nuclei at higher energies, in regions of high 
density, we need an interval of energy /\E centered at E like E » AE >> 1. As 
can be seen from the numerical results, see the top right and the bottom of Figure 1, 
when E » AE, the nearest-neighbor spacing distribution of the vacuum spectrum 
are consistent with a kind of mixing of GOE distribution and GUE, with the same 
characteristic features of random matrix, i.e. level repulsion at short distances and 
suppression of fluctuations at large distances. The histograms at the bottom of the 
Figure 1 gives 2000 consecutive levels centered at differents E values in order to 
show the spectral rigidity, i.e. a very small fluctuation around its average. 

5 A model for the quantum vacuum 

Following [3] we present a model that approximates the above quantum vacuum. 
It could be represented by the non-relativistic s-wave scattering by a spherically 
symmetric potential barrier, V{r) = Vq for r < R and V{r) = for r > R. From 
the stationary scattering solution one obtains the Jost function f+{k) and after in- 
troducing dimensionless parameters, j3 = kR and v = VqR?, its zeros are given by 
the solutions of the complex transcendental equation \J fP' — v cot \/ fi"^ — v = 
for each value of potential strength v. The displacements of the roots, in the 
/3-plane with the variation of the potential strength was shown many years ago by 
Nussensveig |15j. All these zeros are in the lower half of complex /3-plane and lo- 
cated symmetrically about the imaginary axis. The energy/width ratios {AnEn/Gn} 
where calculated (see jSj, Fig. la ) from the Jost zeros for v = 2, in order of dis- 
tance from the origin and show an approximate agreement from the beginning with 
{tn}, for n up to 6 X lO'^. The same Jost zeros for v = 2, after transformation 
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gives the dimensionless widths {G„/47r} in agreement with the global behavior 
of the sequence of primes, from the beginning in the same range (see iSj, Fig. lb ). 
The local behavior, defined as deviations from the average density of the zeta zeros, 
are not obtained by this potential 

The qualitative agreement of the above model suggests the investigation of the 
many-body problem described by one particle being scattered by an effective cutoff 
potential where, by analogy to the mean field, a kind of residual interaction is 
introduced phenomenologically, in order to describe the "vacuum" fluctuation with 
the same statistical distribution of the vacuum energy levels {pntn}- In this case, 
the universality of level fluctuation law of the spectra of different quantum systems 
(nuclei, atoms and molecules) 012] could be understood by the "vacuum" role as a 
dissipative system [Tsj . 

6 Conclusion 

In summary, the spectral interpretation of the Riemann zeta function is associated 
to the vacuum spectrum by means of a variant of the Hilbert-Polya conjecture. The 
distribution of the virtual resonances would reflect a chaotic nature of the quantum 
vacuum. The energy/width ratios of the virtual resonances are associated with 
the nontrivial zeta zeros while the corresponding widths are related to the prime 
sequence. Finally, a weak repulsive cutoff potential is proposed as an approximate 
model for the quantum vacuum. 
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